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We examine the relation between the damping rate of a massless, chiral 
fermion that propagates in a medium, and the rate V of approach to equilib- 
rium. It is proven that these quantities are equal, by showing that they are 
given by the same formula in terms of the imaginary part of the self-energy 
evaluated at the energy of the propagating fermion mode. This result is valid 
provided V is defined by using the appropriate wave functions of the mode. 
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It is well known that a massless fermion that propagates through a medium acquires 
a dispersion relation that differs from that in the vacuum In general, the dispersion 
relation of a fermion with momentum = (e, P) is not given by Ep = P and, in particular, 
Bp is not zero at zero momentum. Writing 

1 

£p =£ r -i (1) 

with e r and 7 being real, the quantity M = s r (P = 0) can be interpreted as an effective 
mass, and 7 as the fermion damping rate f2|. In the framework of Finite Temperature Field 
Theory (FTFT) HQ , the quantity of fundamental interest is the self-energy, from which the 
dispersion relation is determined. 

Besides the fermion damping rate, a distinct concept, also related to the imaginary 
part of the self-energy, has been used in the literature. This quantity, to which we will 
refer as the total reaction rate, is denoted by T and its inverse gives the time scale for a 
distribution function to approach its equilibrium value In Ref. 0, Weldon expressed 

T as a combination of probabilities amplitudes for various processess, weighted by certain 
statistical factors that account for the absortion and emission of the particles by the medium. 
By analysing the one-loop contributions to boson and fermion self-energies, he showed that 
T could be calculated from the imaginary part of the self-energy f7], although no attempt 
was made to connect it with the damping rate. 

In this manner one ends up with two apparently different physical quantities, the damping 
7 and the total reaction rate T, both related to the imaginary part of the self-energy. By 
studying the linear response of a medium to an external field, the equality between T and 
7 was established for the collective color modes in a QCD plasma ||. However, to our 
knowledge the same relation has not been derived for the fermion case, and this fact has 
been the source of some confussion. In a recent article on the subject [||], two different 
formulas are used by the authors to compute the damping and the total reaction rate in 
terms of the self-energy (Eqs. (2.2) and (2.27) of Ref. 0). Our aim in this work is to show 
that no such dichotomy is necessary and, when correctly calculated in terms of the fermion 
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self-energy, the damping and the total rate are given by the same formula. In this fashion, 
the result 7 = T is proved to hold also for fermions. 

Calculations of 7 carried out during the last few years within the realm of hot gauge 
theories, have produced contradictory results, entangled by questions of gauge invariance 
and other problems [|TI]] . It has been pointed out that such ambiguities are a sympton of an 
incomplete calculation, and a resummation of higher order corrections is necessary to get 



the correct result To avoid these complications and make our presentation as clear as 
possible, we examine here the simpler but still instructive situation of fermions interacting 
with a scalar field through a Yukawa interaction. 

We begin by obtaining a formula for 7 in terms of the self-energy, which is valid when 
7 <C e r . This is the physically interesting regime since otherwise the system would be 
overdamped and the concept of a propagating mode is not meaningful. That formula, 
derived by expanding the equation for the dispersion relation up to terms that are at most 
linear in 7 and the absortive part of the self-energy, is the proper expression for the fermion 
damping rate and it corrects or complements formulas that have been quoted in previous 
works PJT^]. We then proceed to demonstrate that the correct formula for 7 coincides with 
T, provided the latter quantity is defined in terms of the transition probabilities calculated 
with the properly normalized spinors that satisfy the effective Dirac equation in the medium 
instead of the free-particle spinors. The equality between the damping and the total rate 
is first illustrated with a one-loop calculation, and then a general proof is given using the 
spectral representation of the fermion self-energy. 

The properties of a particle that propagates through a medium are determined by the 
linear part of the effective (classical) field equation. For a massless fermion with momentum 
p M , the equation, in momentum space, is 

- S e//)^ = , (2) 

with the self-energy £ e // having the form 

S e// = aj> + bji, (3) 



where a and b are functions of the variables 



e = p ■ u , 



P=^e*-p\ (4) 

which give the energy and the magnitude of the 3-momentum P in the rest frame of the 
medium. We have introduced the vector w M representing the velocity 4-vector of the medium 
which in its own rest frame has components (1, 0). 
The self-energy can be decomposed as 

^e// = S r + iSj , (5) 

in terms of the absorptive and dispersive parts 

1 _ 

with S e jj = 7°Egjj7°. Within the framework of the real-time formulation of FTFT, the 
dispersive part of the self-energy is given by 

E r = S llr , (7) 

while the absorptive part can be calculated by any of the following formulas 

v. _ ^ m 

1 — 2np(x) 

= Sl2 (8) 
2in F (x) ' 1 ; 

The E a (, are the elements of the 2x2 self-energy matrix to be computed using the Feynman 
rules of the theory In Eq. (|8|), 

n F (x) = — !— - (9) 
y ' e x + 1 K ' 
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is the Fermi-Dirac distribution for the fermion that propagates in the medium, written in 
terms of the variable 



x = P(p-u-fi), (10) 

where f3 is the inverse temperature and /x is the chemical potential of the fermion. 

The dispersion relations of the particle and antiparticle are determined by requiring 
Eq. (fj) to have non-trivial solutions. The corresponding condition determines also the poles 
of the fermion propagator and can be written in the form 

f(e,P)J(e,P) = 0, (11) 

where 

f{e,P) = (l-a){e-P)-b, 

J(e,P) = (l-a)(e + P)-b, (12) 

with a and b defined in Eq. (|3|). Then, the dispersion relations Bp are obtained as the 
solutions of 

/(e P ,P) = (13) 

and 

7(e P ,P) = 0. (14) 

In general, the solutions Bp are complex, and a consistent interpretation in terms of the dis- 
persion relation for a mode propagating through a medium is possible only if the imaginary 
part of Bp is small compared to its real part. In this case the mode can be visualized as 
a particle with an energy and a damping rate given by the real and imaginary part of Ep, 
respectively |TB . 



Under such assumptions, each one of Eqs. ( [13] ) and ( |HD yields two distinct solutions, 
one with positive energy and the other with negative energy, whose physical interpretation 



has been discussed in detail in Ref. ( [0]). Here we will concentrate on the solution of 
Eq. having a positive real part, which corresponds to the particle mode with energy e r , 
but similar considerations and results apply to the other solutions. 

It is convenient to separate the function / into its dispersive and absorptive parts, and 
write |T3| 



f(e,P) = f r (e,P)+if i (e,P) 



(15) 



with a similar decomposition for a and b. Then 



f r = {l-a r )(e-P)-b r 
fi = -ai{e - P) - hi. 



(16) 



Writing ep in the form of Eq. ([]]), the condition ([H|) becomes 



f r ( £r -il,P)+if i ( £r -il,P) = 



(17) 



which we solve by expanding in powers of 7 and retaining only terms that are at most linear 
in 7 and fi. Thus, e r is determined as the solution of 



/ r (e r ,P) = 0. 



(18) 



while 7 is given by the formula 



where 



7 = fi(e r ,P) 
2 N P 



(19) 



dfr 

de 



(20) 



It is important to observe that the quantity Np coincides with the normalization factor that 
has to be taken into account in order to construct the spinors representing the one-particle 
states, and which need to be included to correctly calculate the probability amplitudes for 
the various processes involving the fermion |I6| . We also want to remark that our Eqs. 
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(|18|) and (|T9|), with the factor Np included, are analogous to those used in the context of 
many-body physics |T7 |. 

The one-particle states have associated the wave functions ul.r satisfying the Dirac 
equation 

- E r ) u L , R = , (21) 

obtained by discarding the absorptive part of the self-energy. While the explicit solutions of 
this equation are more easily worked out in the rest frame of the medium, for our present 
purposes the knowledge of the spinor projection operator will be sufficient. For the negative- 
helicity solution it is given by 

u L u L = 4r L fi , (22) 
i V p 

where, as usual, L = (1 — 7s)/2, Np has been defined in Eq. ( ^0|) and 

n» = i (p" - (e r - P)u») (23) 

is a null vector, whose components in the rest frame of the medium are n M = (1, There 
is a relation for the positive- helicity solution up similar to Eq. ( 122|) . For definiteness, in 
the rest of the article we will work with the negative-helicity solution, but exactly the same 
considerations and results apply to the other one as well. Using 

= a4 + hi , (24) 

and Eq. (p2|), it is easy to verify by direct computation that the formula for the damping 
rate given in Eq. (|19D can be written in the form 

7 = "^Tr m(e r ,P)] 
= -—u L Zi(e r ,P)u L . (25) 

In Eq. ( PS] ) we have indicated explicitly that Sj must be evaluated at e = e r . This a 
convenient formula for 7 that will allow us to identify it with the rate T computed by squaring 



amplitudes and integrating over the phase space with the appropriate statistical weights. 
One should become aware of the general nature of Eq. (|25|) , which we have derived under the 
assumption that the damping is small, but otherwise not tied in anyway to a perturbative, 
or for that matter any other, particular method of calculation of the self-energy. 

If e r and Np are approximated by their value in vacuum (e r ~ P and Np ~ 1), then it is 
easy to see from Eq. (|l^) that 7 ~ — 2Im6(P, P), which coincides with the formula given in 
Ref. ( [T^]) for the damping rate in the high momentum limit (P > T). However, in other 



momentum regimes 

7 ^ -2Im6(e r ,P), 

and the damping rate must be calculated as given in Eq. fllTf ) (or Eq. fl25|)), and by using the 
true dispersion relation in the medium e r . In particular, the formula for the damping rate 
quoted in Eq. (2.27) of Ref. ( 0), which is equivalent to put 7 = — 2lmb(e r , P), does not 
correspond to the imaginary part of the energy momentum-relation and does not coincide 
with the total reaction rate. 

The next step is to establish the relation between Ej and the total reaction rate T, which 
we will do first by using a 1-loop example calculation of Ej. As an instructive example, that 
is free of the ambiguities associated with gauge invariance pointed out in the beginning, we 
adopt the model of a chiral fermion field fi interacting with a scalar field and a massive 
fermion ip according to 



v int 



\^ R f L (f> + h.c. (26) 



The simplest way to proceed with the calculation is to compute E12 and then use Eq. ([8]). The 
relevant diagram is shown in Fig. 1 and, according to the Feynman rules, its contribution 
is given by 

-zE 12 (p) = (zA)(-zA) J -0-iA 21 (p'-p)RiS 12 (p')L. (27) 

The components of the propagator matrices for the massive fermion and the scalar are 
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S 12 (p') = 2m5(p' 2 - m%) [rj F {p') - 6{-p' ■ u)) (/ + m^) 
A 21 (k) = -1-nibik 2 - ml) [r] B {k) + 9{k ■ «)] , 



(28) 



where, in terms of the variables 



x' = (3(p' ■ u- m) , 



(29) 



we have 



Vf(p') — 9(p' • u)n F {x') + 9(— p' ■ u)np{—x') . 
rj B {k) = 9(k ■ u)n B (x lf> ) + 9(-k ■ v^n^-x^) 



(30) 



with 9 being the step function. The fermion distribution rip is given in Eq. ([|), while for 
the scalar 



nsix^) 



e 2 ^ — 1 



(31) 



Momentum as well as charge and lepton number conservation imply that 

X X — X ' (p ^ 

where x has been defined in Eq. (|10|). 

The propagators Si 2 and A 2 i can be rewritten using the identities 

rj F {p' ■ u) — 9(—p' ■ u) = e(p' ■ u)n F (x') , 
rj B (k ■ u) + 9(k ■ u) = e(k ■ u)n B {x (t) )e x ' t ' , 

where e(x) = 9{x) — 9(—x). Substituting the resulting formulas into Eq. (^ 
Eq. (H), we obtain 

A 2 



(32) 



(33) 



and then using 



d p'5{p - m^)5 ip' -p) -m^ 



x e(p' ■ u)e ((p' -p)-u) {n F (x ) + n B {x (j) ))R{j/)' + m^)L. 



(34) 



To derive this result we used the equality 

e^n F (x')n B {x r j>) = n F (x)(n F (x') + n B {x^)) , 

which is a consequence of Eq. fl32|). The integral in Eq. fl34|) can be evaluated for different 
cases, but here we do not go any further in that direction. For our purposes it is more 
convenient to express Eq. (PIP as 

+ 5 (4) (p — k + p') (rfy + Tfy) (jf - m f ) 
+ <5 (4) - A; -p')(l - + n^){j// + m^) 

+ 5 (4) (p + k + p')(l - % + n^W ~ m *l>) ] A ( 35 ) 

where n^, and are the particle density distributions 

1 

U ^ ~ e P{ P '-u- N ,) + l ' 

n * = e /3(fc-u-^) _ 1 ' ( 36 ^ 

and are the respective antiparticle distributions, obtained from by changing the 
sign of the chemical potential. In addition, 

p'» = (E',p'), 
P = (uj k ,k) , 

with 

E' = sjpi* + m\ , 
u k = \jk 2 + ml . 

Now, by means of the projection operators 

u^(p', s)u^{p', s) =ft + m^, 

s 

^ ity(p', s)v f (p', s) =p" -m f , (37) 

s 
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Eq. ( p5[) leads to the relation 

u L Ei{e r ,P)u L = -e r T, (38) 

where we have defined 

F = — f d3p> (Pit) 4 

~ 2e r J (2ir) 3 2u k (2n) 3 2E' [ ' 

x { s^( P + k- p')K(i + + M 1 - \ M i\ 2 

s 

+ 5 {A \p-k + p)[n i ,(l + n^+n 4> (l-n 1 p)}Y,\M n \ 2 

s 

+ 5&(p -k- p')[(l - n^)(l + n^) + |M 7 | 2 

s 

+ 5 (4) (p + A; + p')[(l - n^)(l + n^) + TfyftyE |M /7 | 2 } (39) 



and 



M/ = Xu^{p', s)u L (p) , 

M n = \v4p',s)u l ( P ) . (40) 



The formula for T given in Eq. ( p9| ) is immediately recognized as the total rate for a 
particle of energy e r and momentum P (as seen from the rest frame of the medium) with 
integrations over the phase space weighted by the statistical factors appropriate for each 
process. Mi is the amplitude for fcf) —> i)j or the decay / — > tpcj), while Mjj is the amplitude 
for ftp — >• or f<pip — > 0. The amplitudes for the inverse reactions are given by the complex 
conjugates of Mi and Mn. For specific values of e r and P, some of these processes will 



be kinematically forbidden and will not contribute to T. As is well known |Tj|, for Fermi 
systems the inverse reactions are inhibited as a consequence of the Pauli blocking effect, and 
they contribute additively to the depletion of the state. Therefore, in this case T is given 
by the sum of the rates for the direct and inverse processes, instead of their difference as 
in the bosonic case. Although our expression for F in Eq. ( p9|) has the same form as the 
one used in Ref. [Q, it is important to keep in mind that for us, the propagating fermion is 
represented in the amplitudes Mi u by the spinor ul that obeys the effective Dirac equation 
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in the medium, and not by a free particle spinor. Comparing Eqs. ( p5|) and (|38|) we finally 
obtain the relation 

7 = r. (4i) 

We next show that the relation between the matrix element of the absorptive part Sj 
of the self-energy and the total rate, that was established above by looking at the 1-loop 
calculation, is a general result. For this purpose we recall that, in the real-time formulation 
of FTFT, the self-energy of a fermion in a thermal background is a 2 x 2 matrix whose 
elements are defined by 

z£ 2 i0 - V)cci3 = -(Va(z)rj$(y)) » 

z£l 2 (z - y)a/3 = (Vp(y)Va(z)) , 

-E u (« -y) = E 21 (z - y)6(z° - y°) + E 12 (z - y)6( y ° - z°) , 

-S 22 (^ -y) = E 21 (z - y)6(y° - z°) + E 12 (z - y)6(z° - y°) , (42) 

where t] and rj are the fermion source fields, in terms of which the interaction Lagrangian is 

hnt -IlV + vIl. (43) 



For the scalar model of Eq. (^) rj = X^i/jr. The angle brackets in Eq. ([42]) stand for the 
statistical average which, for any operator O, is defined by 

(0> = iE<W>, (44) 

where 

p = e -P H +Y.A«AQA (45) 

and 

Z = J2{i\p\i), (46) 
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with H being the Hamiltonian of the system. The quantities Qa are the (conserved) charges 
that commute with H, and the a a are the chemical potentials that characterize the compo- 
sition of the background. 

Using the convention of Eq. (fEJ), the amplitude for the decay n + f(p) — > m, where n 
and m label any two states of the system, is 

A ={2tx) a 5^\p + q n -q m ){m\^)u L {p)\n), (47) 

where ul{p) is the properly normalized spinor of the propagating mode /. From this follows 
that the total decay rate, averaged over the initial states of the system, is 

r ^ = \h E I Hr](0)Mp)\n)\ 2 (2^\p + q n - q m )Z n , (48) 

^ ^P n.m 



where Z is defined in Eq. (46) and 



Z n = (n\p\n) 

Here a n is the eigenvalue of the operator a = Y.a°-aQa corresponding to the state |n); i.e., 
a\n) = a n \n). In similar fashion, the total probability for the inverse decay m — > n + f(p) 



is given by 



where 



r ' = 4 A E I (n\u L (pM0)H\ 2 W^ ] (P + Qn- q m )Z m . (51) 

^ ^P n,m 

Since 77 and / have the same quantum numbers, then 

(n\uL(p)7)(0)\m) 7^ 

only for those states such that a m = a n + aj . This relation, combined with the momentum 
conservation condition implied by the delta function in the formulas for Tdj above, allows 
us to replace in Eq. (|5l|) 
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Z m = e~^ u+a fZ ni (52) 

which gives 

T/ = e- x T D , (53) 

or equivalent ly 

(l-n F (x))T I = n F (x)r D , (54) 

where n F and x are given in Eqs. (|9|) and (|10l), with /i = 

In order to establish the relation between Ej and the total rate, we start from the defining 
Eq. ((42j). Inserting a complete set of states J2 \ n ){ n \ between the fields rj(y) and rj(z), we 
have 

iZ 12 (z - y) af3 = \ E e"*^-^-^) x 

(m|^ /3 (0)|n}(n|7/ Q ,(0)|m)Z m , (55) 

from which we immediately obtain 

iu L (p)ll 12 (p)u L (p) = 2p Q T I . (56) 

In similar fashion, 

zE 21 (z - y) aP = ^ E e^"^^) x 

(n|^ a (0)|m)(m|^ / 9(0)|n)2 n , (57) 

so that 

iu L {p)Y, 2l {p)u L {p) = -2p°T D . (58) 
Using these results in the relation Ej = |(E 2 i — Si 2 ) [Eq. (||)] we finally obtain fl9| 

u L {p]Ei{p)u L {p) = -p°T , (59) 

where 
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r = r, + r D . (60) 



On the other hand, recall from Eq. fl2q ) that the damping rate, determined as the imaginary 
part of the dispersion relation, is given by 

7 = -j-u L Y,i(e r ,P)u L . 

where e r is the (real part of the) true dispersion relation of the propagating mode, and ul is 
the corresponding spinor wave function. With the spinor ul thus chosen and with p° = e r , 
T in Eq. fl60|) has the interpretation of the total rate for all the processes involving the 



propagating fermion, so that Eqs. (|59l) and (|25|) establish the equality 7 = T with generality. 

In conclussion, the damping rate, determined from the imaginary part of the dispersion 
relation, coincides with the total rate, provided the latter is calculated with the correct 
wave function of the propagating fermion mode. This result provides a clear and consistent 
interpretation of the absorptive part of the self-energy in terms of a single physical quantity. 
Further, since the damping rate is expressed in terms of transition probabilities, our work 
could serve to consider the problem of the gauge dependence of the fermion damping rate 
from a different point of view. 
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Figure Captions 

Fig. 1. Self-energy diagram for the fermion / due to the interaction \ip R f L 4>. 
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